The Strange Man in Random Networks of Automata 
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We have performed computer simulations of Kauffman's automata on several graphs such as the 
regular square lattice and invasion percolation clusters in order to investigate phase transitions, 
radial distributions of the mean total damage (dynamical exponent z) and propagation speeds of 
the damage when one adds a damaging agent, nicknamed "strange man". Despite the increase in 
the damaging efficiency, we have not observed any appreciable change at the transition threshold 
to chaos neither for the short-range nor for the small-world case on the square lattices when the 
strange man is added in comparison to when small initial damages are inserted in the system. 
The propagation speed of the damage cloud until touching the border of the system in both cases 
obeys a power law with a critical exponent a that strongly depends on the lattice. Particularly, 
we have ckecked the damage spreading when some connections are removed on the square 
lattice and when one considers special invasion percolation clusters (high boundary-saturation clus- 
ters). It is seen that the propagation speed in these systems is quite sensible to the degree of dilution. 
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I. INTRODUCTION 



The study of cellular automata has interested physi- 
cists because of their many applications [l|, Q and their 
contributions to theoretical physics f3(. In fact, problems 
of cellular automata are very related to classical mod- 
els of statistical mechanics (percolation, spin-glass, Ising 
model) . Furthermore cellular automata can be seen as an 
alternative to equations of motion [3], although, nobody 
has yet found a clear way to switch from one description 
to the other. A variety of striking phenomena as stabil- 
ity, criticality, fractals and chaos, can be discussed within 
the technique of automata. In particular, we are inter- 
ested in studying how small failures (damages) produced 
at complex structures of automata propagate through- 
out the system. A biological interpretation can be given 
to the problem of mutations on biological genes. A mu- 
tation could be a change at any particular point (site) 
of the genetic material which affects the time develop- 
ment of the system. The diference between the genet- 
ical material with mutation and without occurrence of 
mutation is described by the Hamming distance of the 
system. Other authors [1, H, 0, have made computer 
simulations finding a phase transition as well as fractal 
behaviour in various structures by inserting a small ini- 
tial damage and then allowing the system to evolve in 
time. But what does happen if one introduces an agent 
that damages the system all the time? We will call this 
agent " strange man" because his rules are randomly up- 
dated during the whole process, i.e., he is a probabilistic 
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automaton among deterministic automata. We gave him 
this nickname in analogy to a social model where each in- 
dividual follows the rules very well, except one. We have 
performed numerical simulations in order to investigate 
phase transitions, radial distributions of the mean to- 
tal damage and propagation speeds of the damage under 
such conditions. We have dealt with well-known networks 
of automata such as regular square lattices (short-range 
case) , invasion percolation clusters M and also the small- 
world topology (small- world case)|10j by making use of 
one of the most important and physically meaningful au- 
tomata model, namely the Kauffman model. The Kauff- 
man model or more generally random Boolean networks 
have got renewed interest because of their recent appli- 
cations to sychronization [llj, stability [13 ■ control of 
chaos and scale-free topologies [T^.ll5|. In our simu- 
lations, despite of the increase in the damaging efficiency, 
we have not observed any appreciable change of the tran- 
sition threshold to chaos neither for the short-range nor 
for the small-world case on the square lattice when the 
strange man is added in comparison to when only small 
initial damages are inserted in the system. Nevertheless, 
the propagation speed of the damage cloud until touch- 
ing the system boundaries in both cases is substantially 
increased and obeys a power-law. For the short-range 
case, we have used an extrapolation technique in order 
to reduce finite-size effects on the square lattice. Further- 
more, we have also checked the damage spreading when 
some connections are removed from the square lattice 
and when one considers invasion percolation clusters with 
high boundary-saturation. It is seen that the propaga- 
tion speed in these systems is quite sensitive with respect 
to dilution. In particular, the damage remains confined 
around the strange man when only few connections are 
removed (typically less than 20%). That shows the ex- 
istence of unstable zones (Tg} in the system, i.e., zones 
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formed by susceptible sites over which the damage can 
spread. 



II. THE KAUFFMAN MODEL 

In the past, KaufFman [l^ introduced networks of 
Boolean automata in order to study the behaviour of 
generic regulatory systems. The basic idea of the Kauff- 
man model is to consider a mixture of all possible binary 
cellular automata. The Kauffman model can be realized 
on a lattice, where Boolean rules are chosen separately 
for each site. Each of N lattice sites hosts a Boolean 
variable cr^ ( spin up or down) which is either zero or 
unity. The time evolution of this model is determined 
by N functions fi (rules) which are randomly chosen for 
each site independently, and by the choice of K input 
sites {Jif(j)} for each site i. Thus the value ai at site i 
for time t+1 is given by: 

a,{t+l) = Ma,,{t),...,aj,At)) i = l,2,...,N. (1) 

Each Boolean function fi is specified, once its value is 
given for each of the 2^ possible neighbour configura- 
tions. A variable ai is said relevant for the spreading 
damage process if it is unstable and/or the state of other 
variables {uj} depend on ai. Moreover if one imposes 
that both the inputs and the chosen Boolean functions 
do not change with the time, we have the quenched Kauff- 
man model. On the other hand, if one admits that both 
change with time, we have the annealed KaufFman model. 
A big difference between the two cases is that in the 
quenched case there are limit cycles and in the annealed 
case not. As the time development in this case is to- 
tally deterministic, and since N different Boolean vari- 
ables can produce 2^ different lattice configurations, we 
must return after at most 2^ time-steps to the previous 
initial configuration. Then the system will repeat the 
same configurations, without ever leaving this limit cycle. 
For the nearest-neighbour Kauffman model on the square 
lattice, the number of relevant limit cycles increases ex- 
ponentially with system size in the non-chaotic phase, as 
shown by Derrida and Flyvbjerg p^ . Kauffman iden- 
tified these different limit cycles with the different cell 
types in our body and found that their number grows as 
^/N for TV interacting genes. Thus the nearest-neighbour 
square lattice is not the most realistic model for such bi- 
ological aspects. The annealed case can be solved ana- 
lytically, whereas for the quenched case only computer 
simulations were performed up to now. 



III. COMPUTATIONAL PROCEDURE 

A standard fashion to implement the Kauffman model 
on computers is introducing a parameter p such that for 
each site on the lattice we select among the 2^^ rules one 
which for each outcome will have spin up with probability 



p. In a computer simulation, first one goes through all N 
sites of the system, and for each site one goes through all 
2^ neighbour configurations, and for each such config- 
uration one determines by drawing one random number 
if its spin will be up or down; if the random number is 
smaller than p then its spin will be up, otherwise it will 
be down. Once one has gone through all neighbour con- 
figurations of that site, then one has fixed the rule for 
that site, and one can go to the next site. After that 
one selects an initial configuration of the Boolean vari- 
ables by randomly assigning to each lattice site a spin 
up or down with equal probability. We will consider 
two systems (replicas), identical in the connections and 
rules, and also identical in the initial configuration of the 
Boolean variables, except that on one of them we put the 
strange man on a central site of the lattice. The number 
of spins which at time t are different between the two 
replicas is called the Hamming distance d(t) or simply 
the damage. For two lattice configurations {ai{t)} and 
{pi{t)}, we have 

m ^ ^T.\''^{t) ~ Mt)i (2) 

i 

and we can define an order parameter ip for the system 
taking in Eq. [2] the limit t ~f oo , that is 

ih = lim d(oo). (3) 

d(0)^0 

Computationally, convergence is typically reached after 
some few thousand time steps. In this way we can study 
both the phase transition and the propagation speed of 
the damage cloud on all networks treated here, by varying 
the value p. Obviously, p and 1—p are statistically equiva- 
lent, so that we do not consider p > 0.50. In order to take 
into account the small-world case, we have also intro- 
duced a rewiring probability q which rewires the K = A 
inputs of each site on the square lattice to any other site, 
not necessarily being its nearest-neighbour. For more 
details about this rewiring procedure as well as a phase 
transiton and fractality in the Kauffman model on small- 
world topology, the reader should consult Ref.pL9]. In 
the present work, we have regarded both the quenched 
and the annealed version. It is worth noting that for the 
annealed version in the infinite range limit, Eq. [5] can be 
iterated [l^l as 

d{t + l) = 2p{l-p)(l-{l~d{t))^), (4) 

where (1 — d{t))^ is the fraction of sites that at the next 
time-step have the same K inputs. 

Invasion percolation clusters were created starting 
from a square lattice and assigning to each lattice site 
a random number between and 1. The algorithm to 
produce such a cluster is as follows: First one chooses 
a central site on the lattice to be the initial seed of the 
cluster. Then one chooses the site with the smallest ran- 
dom number among its nearest-neighbours and adds this 
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FIG. 1: Time evolution of the Hamming distance on a 80x80 ^1^. 2; Time evolution of the Hamming distance on a 80x80 
lattice for different values of w at p = 0.30. lattice for different values of m at p = 0.40. 



site to the cluster. That defines a new set of nearest- 
neighbours for the cluster, i.e., a new boundary of the 
cluster. Then again one connects to the cluster the site 
with the smallest random number belonging to this new 
boundary and so on. The growth process continues un- 
til the cluster touches the borders of the square lattice 
a certain number of times. This multiple touching on 
the lattice borders has the effect of increasing the num- 
ber of connections in the cluster inside the square lattice. 
Such a procedure was necessary to guarantee the damage 
spreading throughout the system. We called the clus- 
ter built with this procedure " high boundary-saturation 
cluster" (HBSC). In our algorithm, the number of touch- 
ing sites on the lattice borders is controlled by the vari- 
able ^ (boundary-saturation). If ^ = 0.80, then 80% of 
the sites of the border of square lattice will be part of 
this HBSC. Once the HBSC has been created, the strange 
man is placed at the initial seed of one of the cluster repli- 
cas and the comparision between replicas is then carried 
out. 

Furthermore we have considered a diluted square lat- 
tice Kauffman model. The inputs of nearest-neighbours 
now vary from K — 1 to K — 4 for each site. That was 
achieved by introducing a connection probability pimk 
for each of the four inputs of the sites and imposing at 
least one input per site. The sites i are treated one by 
one. One goes clockwise drawing a random number for 
each one of its four inputs {(Jj}; if the random number 
is smaller than punk then that input aj is a variable of 
the Boolean function fi [Eq. [1], otherwise it is not. This 
looping procedure is repeated until at least one input is 
a variable of the Boolean function. In this way discon- 
nected sites of the lattice are avoided in our simulations. 

We also define another probability called the strange 
man's behaviour pattern m controlling how many times 
the strange man damages the system per time unit. It 
could be understood as its activity rate. So if one regards 
a strange man with w = 0.50 then for each time step he 
has a 50% chance to randomly change his rule. 



The quenched version was studied both on the HBSC 
and on the square lattice. We have used the annealed 
version only on invasion percolation clusters to calculate 
both the propagation speed of the damage cloud and the 
radial distributions of the mean total damage on these 
structures. In this case the Boolean functions change 
after each time step. 



IV. NUMERICAL RESULTS 

As a preliminary test for about how the strange man 
acts in an automata system, we simulate two replicas of 
a square lattice in order to analyze the temporal evo- 
lution of the Hamming distance up to 1000 time-steps 
for strange men with different behaviour patterns w. In 
Figs[T]and[2l we show the Hamming distance on a square 
lattice with 6400 sites and five different values of at 
p — 0.30 and p ~ 0.40 respectively, averaged over 100 
different initial configurations. As we can see from Figs. 
[1] and [21 changes in the values of w are only relevant in 
the transient state. Even for a rather small value for zu 
as 0.01 the system converges rapidly towards the station- 
ary state. Due to the little impact caused by changing 
137 when one considers the asymptotic time development 
of the system, we will from now on always keep w = l.O 
along this paper. However a variation in w could imply 
a significant change in the touching time on the lattice 
boundaries, since this happens in the transient state. 

Another interesting issue is to search for a phase tran- 
sition on the square lattice when one introduces a strange 
man in the system. This is done for a 500x500 lattice by 
calculating the order parameter ^ averaged over 100 runs 
with up to 1000 time-steps, where we have also studied 
the impact caused when one cuts randomly some of their 
connections, as shown in Fig. [3l Two conclusions can 
be drawn: On one hand, the continuous activity of the 
strange man does not affect the value of the critical point 
which for the square lattice is at 0.31 according to Ref. 
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FIG. 3: Phase transition on a 500x500 lattice for three dif- 
ferent values of the probabilities punk The critical point at 
Piink = 1-0 is 0.31, at piink = 0.90 is 0.37 and at pii„k = 0.80 
one no longer has a phase transition. 
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FIG. 4: Phase transition on a HBSC built on a 100x100 lat- 
tice for three different values of the variable (, in the quenched 
case. Particularly, at ^ = 0.90 we show our result for the an- 
nealed case as well. At ^ = 0.90 one has pc = 0.26 (annealed 
case) and pc = 0.36 (quenched case). 
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FIG. 5: Log-log plot of ^ vs (p — pc) for a HBSC built on 
a 100x100 lattice for both the annealed and the quenched 
version. 



. The same thing also happens if one regards the small- 
world case jlE]. On the other hand, the damage spread- 
ing on square lattices is quite sensitive to the removal 
of connections. That is due to the existence of unsta- 
ble regions on the lattice (unstable cores). As we can see 
from Fig. 3, taking a connection probability punk = 0.90, 
the value of the critical point on the lattice is pushed to 
p — 0.37. Yet more dramatic, at pimk = 0.80 we do not 
find a chaotic phase any more. 

We have done the same analysis for the order param- 
eter V' on a HBSC built on a 100x100 lattice, where N 
in Eq. [5] is the number of sites of a typical cluster. Our 
results are shown in Fig. |4] for three different values of 
^ (see previous section). Particularly, at ^ = 0.90 we 
also present a simulation for the annealed version on this 
cluster. For a boundary-saturation ^ = 0.90, we can see 



a phase transition around Pc — 0.36 in the quenched ver- 
sion and around Pc — 0.26 in the annealed version. In 
Fig. [5] we show a log- log plot of ip vs {p — Pc)- Our 
estimates for the critical exponents (3 for both versions 
are: 

13 ^ Q 42 ± 0.02 (5) 

^quenched ^ Q 47 _^QQ2 (6) 

Taking advantage of the strange man's high efficiency 
to provoke damage, we have performed several numerical 
calculations in order to evaluate the propagation speed 
of the damage v by measuring the time it takes to touch 
the system boundaries both on the square lattice and 
on the HBSC. The propagation speed on square lattices 
of size L X L was calculated for both the short-range 
case {q — 0) and the small-world case {q ^ 0) in their 
quenched version, where for the short-range case we ap- 
plied an extrapolation technique in order to reduce finite- 
size effects. This extrapolation was achieved by analysing 
how the propagation speed of damage v depends on the 
reciprocal of the lattice size (1/i) when one takes the 
limit L —t oo for every value of p considered. Fig. [5] 
shows the propagation speeds for the short-range case 
for L=100, 300 , 700 and L — oo (extrapolated data) 
averaged over 100 different runs for each parameter p, 
while Fig. [7] shows the propagation speeds for the small- 
world case at g = 0.05 for lattice sizes L=100, L=200 and 
L=300 averaged over 1000 different runs. It is somewhat 
surprising to observe in Fig. [7] a considerable increase of 
the propagation speed compared to the case plotted in 
Fig. [B] as well as a clear dependence of speed on the sys- 
tem size. One can observe an almost constant propaga- 
tion speed above p = 0.45. In other words, it makes little 
difference for the system dynamics if one is at p = 0.45 
or at p = 0.50. We have built clusters with boundary- 



5 



0.8 

0.7 

T3 0.6 
(D 
(D 

Q. 0.5 
W 

I 0.4 

1,0.3 
IB 

9- 0.2 

o 

^ 0.1 





8 ^ 



0.3 0.35 0.4 

(P) 



+ « Q 



O L=100 
A L=300 
V L=700 
+ L= oo 



0.45 



0.5 



FIG. 6: Propagation speed of the damage in the short-range 
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an extrapolation of these data {L — > oo). 
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FIG. 8: Propagation speed of the damage on the HBSC at 
5 = 0.90 inside a 100x100 lattice both in the quenched and 
the annealed case. 



FIG. 9: Log-log plot of v vs (p — pc) for the extrapolated data 
from Fig. [S]in the short-range case (g = 0.0). 



saturation ^ = 0.90 inside square lattices of size 100 x 100 
and made simulations over 1000 different HBSCs both 
for the quenched version and the annealed version. Fig. 
[5] shows our results for the propagation speed on HBSC 
at values above and a little below the critical point. The 
major reason why the propagation speed does not vanish 
at all for values of p a little below the critical point is the 
existence of finite-size effects. Moreover it is interesting 
to observe that for the annealed version one has a higher 
propagation speed compared to the quenched version. 

To estimate the critical exponents aq of the speed for 
the square lattice, we have made a log-log plot of the 
data from Figs. [5] and [7] by taking pc = 0.31 as critical 
point for the short-range case {q — 0.0) and Pc — 0.15 as 
critical point for the small- world case {q = 0.05). Fig. [5] 
shows the log-log plot of v vs (p—pc) for the extrapolated 
data from Fig. [SI our estimate for the critical exponent 
of the speed at {q = 0.0) is w 0.27, while in Fig. [TO] a 
collapse of the data from Fig. [7| is shown by plotting the 



log-log plot ofv/L vs (p—pc), our estimate for the critical 
exponent of the speed at {q — 0.05) is aoos ~ 0.62. In 
addition, we also calculated the critical exponent a of the 
speed on the HBSC at ^ = 0.90 through a log-log plot of 
the data from Fig. [5| as shown in Fig. [TT] Our results 
for the critical exponents a of the speed on the HBSC 
are nearly the same for the quenched and the annealed 
version. For the quenched version we obtain a « 0.30 
whereas for the quenched version we obtain a « 0.33 
where we use Pc — 0.36 and pc — 0.26 as their respective 
critical points. 

But how does the damage cloud emitted by the strange 
man evolve with time when one takes into account dif- 
ferent values of p7 To answer this question, we have per- 
formed simulations in order to determine the dynamical 
exponent z. First one determines the radial distribution 
of the mean total damage Mt- This is defined as the 
mean number of times that the sites along a straight line 
starting from the strange man have been flipped dur- 
ing the whole time development after the transient time. 
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FIG. 10: Log-log plot of v/L vs {p — pc) of the data plotted 
in Fig. [3 in the small- world case (q = 0.05). 




FIG. 12; Semi-log plot of Mt vs the distance r from the 
strange man on the square lattice for three different values of 
p above the critical point, namely p= 0.32, 0.35 and 0.45. 
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FIG. 11: Log-log plot of V vs (p — pc) for the data from Fig. 
[S]for both the quenched and the annealed version on HBSC 
at 5 = 0.90 inside a 100x100 lattice . 



is defined by 



(5 cx (p - PcY 



(8) 



where S is the inverse characteristic length and is the 
critical point. One can determine the dynamical expo- 
nent z by a log-log plot of S vs {p ~ Pc) as shown in 
Figs. [13] and [14] for six different values above the crit- 
ical point on the 200x200 lattice and on the HBSC at 
^ = 0.90 inside a 100x100 lattice, respectively. Our esti- 
mate for the dynamical exponent z on the square lattice 
is z w —0.156, whereas on the HBSC we find 



quenched 



^annealed 



-0.165 
-0.188. 



(9) 
(10) 



The negative sign of z above means that damage de- 
creases as one gets farther from the critical point pc- 



That was done by means of a vectorial storage of the 
total damage for each site along a straight line starting 
from the strange man on 200x200 lattices and on HBSCs 
built inside 100x100 lattices over 1000 difi'erent configu- 
rations for different values of p above the critical point. 
On the HBSC both quenched and annealed versions were 
considered with their respective critical points. Moreover 
in order to get rid of the transient regime we iterated the 
systems until about twice their respective touching times 
for each value of p. Fig. [T2| shows the semi-log plot of Mt 
against the distance r from the strange man for three dif- 
ferent values of p above of the critical point on a 200 x 200 
lattice. One can define the inverse characteristic length 
5 for each value of p above the critical point by fitting 
Mt to 

MT{r) = Aexp((5r - </)), (7) 
where A and are constants. The dynamical exponent 



V. SUMMARY 

We have reported some results from simulations to 
study the damage propagation in the Kauffman model 
on different graphs such as square lattices and invasion 
percolation clusters using a damaging agent that changes 
its rule all the time. In this paper, we have looked for 
critical points by analysing the order parameter for each 
case. We have not observed any changes at the criti- 
cal points for both the short-range and small-world case 
on the square lattice when one introduces small initial 
damage on the lattice. However the propagation speed is 
increased due to the strange man's power to damage the 
system. This has decreased our computational effort be- 
cause less time is necessary to calculate the propagation 
speed as well as the radial distribution of the mean total 
damage on these graphs. The critical exponents for the 
propagation speed as well as the dynamical exponents 
have been calculated with reasonable precision for both 
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FIG. 14: Log-log plot of S va {p — Pc) on the HBSC inside a 
FIG. 13: Log-log plot of <5 vs (p - p.) on a 200x200 lattice. loOxlOO lattice for both quenched and annealed versions. 
The slope gives the dynamical exponent z. 



the short-range case on the square lattice and on invasion 
percolation clusters. 

We found for the short-range case on the square lat- 
tice a critical exponent of the speed ao ~ 0.27 and a 
dynamic exponent z ~ —0.156. For the small-world case 
on the square lattice we obtain aoos ~ 0.62. On in- 
vasion percolation clusters we obtain critical exponents 
of the speed a w 0.30 for the quenched version and 
a « 0.33 for the annealed version and dynamical expo- 
nents z q"™ched ^ _o.l65 and ^^nneaied ^ _o.l88. Fur- 
ther we observe a high sensitivity with respect to damage 
propagation when one cuts randomly some connections 
on the square lattice or when one considers invasion per- 
colation clusters with high boundary-saturation. 



In addition it would also be interesting to look how 
the strange man perturbs the basins of attraction in the 
quenched Kauffman model on the graphs studied here. 
This analysis is important to better understand the dy- 
namics of the system. That will be the subject of future 
investigations. 
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